Normal mode splitting in a coupled system of nanomechanical oscillator and 

parametric amplifier cavity 
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We study how an optical parametric amplifier inside the cavity can affect the normal mode 
splitting behavior of the coupled movable mirror and the cavity field. We work in the resolved 
sideband regime. The spectra exhibit a double-peak structure as the parametric gain is increased. 
Moreover, for a fixed parametric gain, the double-peak structure of the spectrum is more pronounced 
with increasing the input laser power. We give results for mode splitting. The widths of the split 
lines are sensitive to parametric gain. 

PACS numbers: 42.50.Wk, 42.50.Pq, 42.50.Lc 



I. INTRODUCTION 

Recently there has been a major effort in applying 
many of the well tested ideas from quantum optics such 
as squeezing, quantum entanglement to optomechanical 
systems which are macroscopic systems. Thus observa- 
tion of entanglement [J, [2, [1,13, S 11] , squeezing 0, [1| etc 
in optomechanical systems would enable one to study 
quantum behavior at macroscopic scale. This of course 
requires cooling such systems to their ground state and 
significant advances have been made in cooling the me- 
chanical mirror to far below the temperature of the en- 
vironment [i, [13, [H [il [H, [13, El. Further it has 
been pointed out that using optical back action one can 
possibly achieve the ground state cooling in the resolved 
sideband regime where the frequency of the mechanical 
mirror is much larger than the cavity decay rate, that is 

uj,n > K iii,[i3,[il. 

Another key idea from quantum optics is the vacuum 
Rabi splitting [l^ [13] which is due to strong interaction 
between the atoms and the cavity mode. The experi- 
mentalists have worked hard over the years to produce 
stronger and str ong er couplings to produce larger and 
larger splittings [U [13, [11] . Application of these ideas 
to macroscopic systems is challenging as well. In a recent 
paper Kippenberg et al. proposed the possibility of 
normal mode splitting in the resolved sideband regime 
using optomechanical oscillators. In this paper, we pro- 
pose placing a type I optical parametric amplifier inside 
the cavity to increase the coupling between the movable 
mirror and the cavity field, and this should make the ob- 
servation of the normal mode splitting of the movable 
mirror and the output field more accessible. 

The paper is structured as follows. In Sec. II we 
present the model, derive the quantum Langevin equa- 
tions, and give the steady-state mean values. In Sec. Ill 
we present solution to the linearized Langevin equations 
and give the spectrum of the movable mirror. In Sec. 
IV we analyse and estimate the amount of the normal 
mode splitting of the spectra. In Sec. V we calculate the 
spectra of the output field. In Sec. VI we discuss the 
mode splitting of the spectra of the movable mirror and 



the output field. 



II. MODEL 



The system under consideration, sketched in Fig. [U is 
an optical parametric amplifier (OPA) placed within a 
Fabry-Perot cavity formed by one fixed partially trans- 
mitting mirror and one movable perfectly reflecting mir- 
ror in equilibrium with its environment at a low tem- 
perature. The movable mirror is treated as a quantum 
mechanical harmonic oscillator with effective mass m, 
frequency Wm, and energy decay rate 7m . An external 
laser enters the cavity through the fixed mirror, then the 
photons in the cavity will exert a radiation pressure force 
on the movable mirror due to momentum transfer. This 
force is proportional to the instantaneous photon number 
in the cavity. 



cavity axis 
M ► 



OPA 



fixed mirror movable mirror 

FIG. 1: Sketch of the studied system. The cavity contains 
a nonlinear crystal which is pumped by a laser (not shown) 
to produce parametric amplification and to change photon 
statistics in the cavity. 

In the adiabatic limit, the frequency ujm of the movable 
mirror is much smaller than the free spectral range of the 
cavity ^ (c is the speed of light in vacuum and L is the 
cavity length), the scattering of photons to other cavity 
modes can be ignored, thus only one cavity mode lUc is 
considered [H, ■ The Hamiltonian for the system in a 
frame rotating at the laser frequency uj^ can be written 
as 



H 



+ine{c^ -c) + ihG{e'''c^^ - e^'^c^). 



(1) 
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Here Q and P are the dimensionless position and mo- 
mentum operators for the movable mirror, defined by 

Q - and P = with [Q,P] = 2z. In 

Eq. IH), the first term is the energy of the cavity field, 
ric = c^c is the number of the photons inside the cav- 
ity, c and c' are the annihilation and creation operators 
for the cavity field satisfying the commutation relation 
[c, c^] = 1. The second term comes from the coupling of 
the movable mirror to the cavity field via radiation pres- 
sure, the dimensionless parameter y = —^\ hr-^ — is 
the optomechanical coupling constant between the cav- 
ity and the movable mirror. The third term corresponds 
the energy of the movable mirror. The fourth term de- 
scribes the coupling between the input laser field and the 
cavity field, s is related to the input laser power p by 

£ — ^ where k is the cavity decay rate. The last 

term is the coupling between the OPA and the cavity 
field, G is the nonlinear gain of the OPA, and 9 is the 
phase of the field driving the OPA. The parameter G is 
proportional to the pump driving the OPA. 

Using the Heisenberg equations of motion and adding 
the corresponding damping and noise terms, we obtain 
the quantum Langevin equations as follows, 

Q — ^mPi 



P = 2uJmX1^c - i^mQ - ImP + 

c = —i{ujc — ojl — ^mXQ)c + £ + 2Ge'^c 



t _ 



= i{bJc - UJL - uJmXQ)c^ + £ + 2Ge - kc^ + ^/2k.c\^ 

(2) 

Here we have introduced the input vacuum noise operator 
dn with zero mean value, which obeys the correlation 
function in the time domain 12711 



{5c,^{t)5cUf))^5{t-t'l 



{dcUt)ScUt')) = {5clJt)Sc,r,{t')) - 0. 



(3) 



The force ^ is the Brownian noise operator resulting from 
the coupling of the movable mirror to the thermal bath, 
whose mean value is zero, and it has the following corre- 
lation function at temperature T [28j : 



TT UJm 



c^e— l-i-coth(-l^) di,, 

(4) 

where ks is the Boltzmann constant and T is the ther- 
mal bath temperature. Following standard methods from 
quantum optics [2^ , we derive the steady-state solution 
to Eq. ^ by setting all the time derivatives in Eq. ^ 
to zero. They are 



Ps ^0,Qs^ 2x\c 



,2 C 

S I ^ 



iA + 2Ge*^ 
+ A2 - 4G2 



where 



LUc-UJL- i^mXQs 



£, (5) 



(6) 



is the effective cavity detuning, depending on Q^. The Qs 
denotes the new equilibrium position of the movable mir- 
ror relative to that without the driving field. Further Cg 
represents the steady-state amplitude of the cavity field. 
From Eq. (O and Eq. ([6]), we can see A satisfies a fifth 
order equation, it can at most have five real solutions. 
Therefore, the movable mirror displays an optical multi- 
stable behavior 31, 32], which is a nonlinear effect 
induced by the radiation-pressure coupling of the mov- 
able mirror to the cavity field. 



III. RADIATION PRESSURE AND QUANTUM 
FLUCTUATIONS 

In order to investigate the normal mode splitting of the 
movable mirror and the output field, we need to calculate 
the fiuctuations of the system. Since the problem is non- 
linear, we assume that the nonlinearity is weak. Thus we 
can focus on the dynamics of small fluctuations around 
the steady state of the system. Each operator of the sys- 
tem can be written as the sum of its steady-state mean 
value and a small fluctuation with zero mean value, 

Q^Qs + SQ, P = P, + SP, c = c, + Sc. (7) 

Inserting Eq. fT]) into Eq. then assuming \cs\ ^ 1, 
the linearized quantum Langevin equations for the fluc- 
tuation operators take the form 

SQ ^ UmSP, 

SP ^ 2uJrnXic*sSc + Cs5c^) - LU„iSQ - JmSP + 

Sc = - {k + iA)Sc + iuj„iXCsSQ + 2Ge'^(5c^ + V2k(5q„, 

(5ct = -(k - iA)Sc^ - iujmXc*JQ + 2Ge-'%c + V2^(5cL. 

(8) 

Introducing the cavity field quadratures Sx — 5c + Sc^ 
and Sy = i{6c^ — 6c), and the input noise quadratures 
Sxin = Scin + (5c-„ and Syin = - (5ci„), Eq. ^ can 

be rewritten in the matrix form 



f{t)^Afit) + 7^it), 



(9) 



in which f{t) is the column vector of the fluctuations, 
ri{t) is the column vector of the noise sources. Their 
transposes are 



(10) 



f{t f ^{SQ,5P,Sx,5y), 
ri{tY = {Q,^,,^/2K5x^n,\/2K5y^n)■, 
and the matrix A is given by 

/ cj„ 

-Um -Irn >^mX{Cs + C* ) -iuJmX{Cs " C* ) 

iujrnX{cs — c*g) 2G COS 6 — K 2Gsin0-|-A 
uimX{cs+c*) 2Gsin6'-A -(2Gcos6' + k) 

(11) 



V 
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The system is stable only if all the eigenvalues of the ma- 
trix A have negative real parts. The stability conditions 
for the system can be derived by applying the Routh- 
Hurwitz criterion [H, [s^l • This gives 



2k{k^ ~ + A2 + 2k7„0 + 7™ (2^7™ + ajfj > 0, 

+n^m{{^' - 4G2 + A2)2 + (2«7„ + 7^) 
X - 4G2 + ^2) ^ ^2^[2(k2 + 4G2 - A2) 
+ 2k7„,]} > 0, 
k2 _ 4G2 + A2 - 4tJ™x'[|c.pA + iG(c2e-*» 



> 0. 

(12) 

All the external parameters must be chosen to satisfy the 
stability conditions (H^)- 

Taking Fourier transform of Eq. ([8]) by us- 

ing fit) = ^ fico)e-^-'dLU and p (t) = 

2^ /_T /n-'^)e-^"*dc., where /t(-c.) = [/(-a;)]t, then 
solving it, we obtain the position fluctuations of the mov- 
able mirror 



+2Ge-'''cs]6c^n{io) + [{k + i{A - Uj))Cs 
+2Ge^'^cl]Scl{~Lo)} 

+ [(^_,^)2+A2_4G2]^M], 



(13) 



where 
d{uj) 



4u;lx'[A\cs\^+^G{cle-^' 



4G2 



In Eq. (|13p . the first term proportional to x originates 
from radiation pressure, while the second term involving 
^(w) is from the thermal noise. So the position fluctua- 
tions of the movable mirror are now determined by radia- 
tion pressure and the thermal noise. In the case of no cou- 
pling with the cavity field, the movable mirror will make 
Brownian motion, 5Q{uj) = lu„i£,{uj) / (lu'^ — uP' — 17^^^), 
whose susceptibility has a Lorentzian shape centered at 
frequency iOm with width 7„j . 

The spectrum of fluctuations in position of the movable 
mirror is defined by 



-((5g(a;)<50(r!)) + {bQ(Sl)8Q[u:))) = 2'^Sq(uj)8{uj + f)). 

(15) 

To calculate the spectrum, we require the correlation 
functions of the noise sources in the frequency domain, 

(<5c„(w)54(-r!)) = 27r(5(Lj + fi), 

{i{uj)m)) = ^^^^ [1 + ^otHw^)] ^{^ + 



Substituting Eq. US]) and Eq. HH) into Eq. ([151), ^e 
obtain the spectrum of fluctuations in position of the 
movable mirror [ssl ] 

+2Ge'''cf{K - iA) + 2Ge-'^cl{K + iA)] 
+2^a;[(A2 + - uj^ - AG^f + Ak'uj''] 
xcoth(^)}. 

(17) 

In Eq. (jl7p . the first term involving x arises from radia- 
tion pressure, while the second term originates from the 
thermal noise. So the spectrum Sq{lo) of the movable 
mirror depends on radiation pressure and the thermal 



IV. NORMAL MODE SPLITTING AND THE 
EIGENVALUES OF THE MATRIX A 

The structure of all the spectra is determined by the 
eigenvalues of iA (Eq. pT|) ) or the complex zeroes of the 
function d{u}) defined by Eq. p^ . Clearly we need the 
eigenvalues of iA as the solution of (Eq. ^) in Fourier 
domain is /(oj) = iiuj — iA)^^ri{Lo). Let us analyse the 
eigenvalues of Eq. pTjl . Note that in the absence of the 
coupling x=0, the eigenvalues of iA are 



±1/ UJ. 



S±VA2 -4G2 



(18) 



Thus the positive frequencies of the normal modes are 

given by VA^ - 4G2, ^uj^-^-^ (A > 2G,^„ > Y)- 
The case that we consider in this paper corresponds to 



> 2^; A > 2G;k > 7„ 



> K. 



(19) 



The coupling between the normal modes would be 
most efficient in the degenerate case i.e. when a;,„ = 
-\/A2 — 4G2. It is known from cavity QED that the 
normal mode splitting leads to symmetric (asymmetric) 
spectra in the degenerate (nondegenerate) case, provided 
that the dampings of the individual modes are much 
smaller than the coupling constant. Thus the mechanical 
oscillator is like the atomic oscillator, the cavity mode 
in the rotating frame acquires the effective frequency 
-y/A^ — 4G2 which is dependent on the parametric cou- 
pling. An estimate of the splitting can be made by using 
the approximations given by Eq. (|19p and the zeroes of 
djoj). We find that the frequency splitting is given by 



^-^±J{^^^^^Y + ^^W. (20) 



(16) 



where we have defined 



c^„x1c.f [A -I- 2Gsin(0 - 2^)], e^'^ = cl/\c 



(21) 
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It should be borne in mind that Cs , A etc are dependent 
on the parametric couphng G. The sphtting is deter- 
mined by the pump power, the couphngs x and G. 

The parameters used are the same as those in the 
recent successful experiment on optomechanical normal 
mode splitting the wave length of the laser A = 

2ttc/ujl = 1064 nm, i = 25 mm, m = 145 ng, k = 
2tt X 215 X 10^ Hz, cj„ = 27r x 947 x 10^ Hz, T = 300 
mK, the mechanical quality factor Q = oJmHm = 6700, 
parametric phase 9 — n/A. And in the high temperature 
limit kgT 3> fiiOm, we have coth{huj/2kBT) sa iksT /huj. 

Figure [2] shows the roots of d{Ld) in the domain 
Re(a;) > for different values of G. Figure [3] shows 
imaginary parts of the roots of d{uj) for different values 
of G. The parametric coupling affects the width of the 
lines and this for certain range of parameters aids in pro- 
ducing well split lines. One root broadens and the other 
root narrows. The root that broadens is the one that 
moves further away from the position for G = Q. 
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FIG. 2: (Color online) The roots of d{Lu) in the domain 
Re(cij) > as a function of parametric gain, p — 6.9 mW 
(dotted line), p = 10.7 mW (dashed line). Parameters: the 
cavity detuning A = LUrn- 
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FIG. 3: (Color online) The imaginary parts of the roots of 
d{uj) as a function of parametric gain, p — 6.9 mW ( dotted 
line), p — 10.7 mW (dashed line). Parameters: the cavity 
detuning A = ujm- 



V. THE SPECTRA OF THE OUTPUT FIELD 

In this section, we would like to calculate the spectra 
of the output field. The fluctuations Sc{uj) of the cavity 
field can be obtained from Eq. ([5]). Further using the 
input-output relation [sll Couti^) = \/2kc{uj) — Cin{uj), 
the fluctuations of the output field are given by 

5cout{Lo) = V{uj)i{Lo) + E{uj)5c,n{uj) + F {lo)6c\J~uj) , 

(22) 

where 



[-'-^^{[^-^{u^ + ^)]cs 
2Ge-*^cJ 



-2Ge*«c*}{[K-i(u;-f A)]c: 



-l-K - i(w -f A)] - 1, 
(.-»c.)4a^-4G^ h^W'^t^ - + 
-2Ge*«c*}{[K - i{oj - A)]c, 2Ge*«c*} 
+2Ge'% 



In Eq. 



(23) 



the first term associated with ^{u) stems 
from the thermal noise of the mechanical oscillator, while 
the other two terms are from the input vacuum noise. So 
the fluctuations of the output field are infiuenced by the 
thermal noise and the input vacuum noise. 
The spectra of the output field are defined as 



{5cl^^{-Vt)5Cout{i^)) = 2TTScouti'^)S{uJ + Q), 



(24) 



{5yout{^)5yout{i-^)) = 2tt Syout{'^)5{uj + Q). 

where 5xout{^) and 6yout{^) are the Fourier transform of 
the fluctuations Sxout{t) and Syout{t) of the output field 
, which are defined by Sxout{t) = Scout{t) + 6cl^f.{t) and 
Syoutit) = «[^cj,„t(t) - Scoutit)] [li]. Here 5'co«t(w) de- 
notes the spectral density of the output field, Sxout{i^) 
means the spectrum of fluctuations in the x quadrature 
of the output field, and Syout{^) is the spectrum of fluc- 
tuations in the y quadrature of the output held. 

Combining Eq. dUl), Eq. (gS)), and Eq. ([Ml), we obtain 
the spectra of the output field 

V*{uj)V{u) X 2^c.hl + coth(5|^)] 



Syout{i^) 



^2kBTl 

+F*{uj)F{uj), 

[V{-uj) + V*{uj)][V{uj) + V*{-cu)] 

x25:^[-i + coth(^)] 

+ [E{~uj) + F* (uj)] [F{uj) + E* {-uo)l 
-[V*{uj)-V{-uj)][V*{-u)-V{u)] 
x25:^[-l + coth(^)] 
~[F*{oj)^E{-u)][E*{-u)^F{u,)]. 



(25) 
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From Eq. (|25|) , it is seen that any spectrum of the out- 
put field includes two terms, the first term is from the 
contribution of the thermal noise of the mechanical os- 
cillator, the second term is from the contribution of the 
input vacuum noise. 

We note that the spectra Sq{uj), Scout{'^), Sxout{'^), 
and Syout{uj) are determined by the detuning A, para- 
metric gain G, parametric phase 9, input laser power p, 
and cavity length L. In the following we will concentrate 
on discussing the dependence of the spectra on paramet- 
ric gain and input laser power. 
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VI. NUMERICAL RESULTS 

In this section, we numerically evaluate the spectra 
Sq^lo), Scoiit{io), Sxout{io), and Syout{u)) given by Eq. 
(fT7)l and Eq. to show the effect of an OPA in the 

cavity on the normal mode splitting of the movable mir- 
ror and the output field. 

We consider the degenerate case A = cj^ for G ~ Q, 
and choose p = 6.9 mW to satisfy the stability condi- 
tions (12), parametric gain must satisfy G < 1.62k. Fig- 
ures |3]-[7]shows the spectra Sq{lli), Scouti^), S^outi^), 
and Syoutito) as a function of the normalized frequency 
uj/ujm for various values of parametric gain. When the 
OPA is absent (G = 0), the spectra barely show the nor- 
mal mode splitting. As parametric gain is increased, the 
normal mode splitting becomes observable. This is due to 
significant changes in the line widths and position. When 
G = 1.3k, two peaks can be found in the spectra. Note 
that the separation between two peaks becomes larger as 
parametric gain increases. The reason is that increasing 
the parametric gain causes a stronger coupling between 
the movable mirror and the cavity field due to an increase 
in the photon number in the cavity. The values of inter- 
cavity photon number |csp are 2.68 x 10^, 4.30 x 10^, 
5.65 X 10^ for G — 0, 1.3k, and 1.45k respectively. We 
have examined the contributions of various terms in Eq. 
(I^SI) to the output spectrum. The dominant contribution 
comes from the mechanical oscillator. Note further the 
similarity j36j of the spectrum of the output quadrature 
y (Fig- El to the spectrum of the mechanical oscillator 
(Fig. 3]). It should be borne in mind that the strong 
asymmetries in the spectra for G ^ arise from the fact 
that by fixing A at a;,„, the frequencies of the cavity mode 
and the mechanical oscillator do not coincide if G 7^ 0; 
X — 0- Besides the damping term k, k being not neg- 
ligible compared to A, also contributes to asymmetries. 

Now we fix parametric gain G = 1.3k, and choose A = 
•y/aj^7+4G^: the input laser power must satisfy p < 55 
mW. The spectrum Sq{uj) as a function of the normalized 
frequency uj/ujm for increasing the input laser power is 
shown in Fig. [51 As we increase the laser power from 0.6 
mW to 10.7 mW, the spectrum exhibits a doublet and 
the peak separation is proportional to the laser power, 
because the coupling between the movable mirror and 



FIG. 4: (Color online) The scaled spectrum Sq{lo)x 7™ versus 
the normalized frequency lo/uJm for different parametric gain. 
G= (solid curve), 1.3k (dotted curve), 1.45k (dashed curve). 
Parameters: the cavity detuning A = u)m, the laser power 
p = 6.9 mW. 




FIG. 5: (Color online) The spectrum Scout{ijj) versus the nor- 
malized frequency ui/uim for different parametric gain. G= 
(solid curve), 1.3k (dotted curve), 1.45k (dashed curve). 
Parameters: the cavity detuning A = Um, the laser power 
p = 6.9 mW. 




FIG. 6: (Color online) The spectrum Sxout{i^) versus the nor- 
malized frequency u;/uJm for different parametric gain. G= 
(solid curve), 1.3k (dotted curve), 1.45k (dashed curve). 
Parameters: the cavity detuning A = Um, the laser power 
p = 6.9 mW. 



6 




FIG. 7: (Color online) The spectrum Syout{i^) versus the nor- 
malized frequency uj/uim for different parametric gain. G= 
(solid curve), 1.3k (dotted curve), 1.45k (dashed curve). 
Parameters: the cavity detuning A = cOm, the laser power 
p — 6.9 mW. 



without OPA {G = 0), the spectrum Sq{uj) as a function 
of the normahzed frequency to/wm for increasing the in- 
put laser power at A = is plotted in Fig. [S) We can 
see if the laser power is increased from 0.6 mW to 10.7 
mW, the spectrum also displays normal mode splitting. 
However the normal mode with OPA (Fig. [5]) are more 
pronounced than that in the absence of OPA (Fig. 
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the cavity field for a given parametric gain G is increased 
with increasing the input laser power due to an increase 
in photon number. 





tJ=0.6 mW 




SJ=6.9 raw 




,J=10.7 mW 



FIG. 8: (Color online) The scaled spectrum Sq{uj) x 7„i ver- 
sus the normalized frequency cj/uim, each curve corresponds 
to a different input laser power. p= 0.6 mW (solid curve, left- 
most vertical scale), 6.9 mW (dotted curve, rightmost vertical 
scale), 10.7 mW (dashed curve, rightmost vertical scale). Pa- 
rameters: the cavity detuning A = \/ lu'^ + AG'^ , parametric 
gain G — 1.3k. 

For comparison, we also consider the case of the cavity 



FIG. 9: (Color online) The scaled spectrum Sq{uj) x 7™ ver- 
sus the normalized frequency uj/um, each curve corresponds 
to a different input laser power, p— 0.6 mW (solid curve, 
leftmost vertical scale), 6.9 mW (dotted curve, rightmost ver- 
tical scale), 10.7 mW (dashed curve, rightmost vertical scale). 
Parameters: the cavity detuning A = ujm, parametric gain 
G = 0. 



VII. CONCLUSIONS 

In conclusion, we have shown how the normal mode 
splitting behavior of the movable mirror and the output 
field is affected by the OPA in the cavity. We work in the 
resolved sideband regime and operate under the stabil- 
ity conditions (12). We find that increasing parametric 
gain can make the spectra S'q(w), Scoutii^), Sxouti'jj), and 
Syout{^) evolve from a single peak to two peaks. Further- 
more, for a given parametric gain, increasing input laser 
power can increase the amount of normal mode splitting 
of the movable mirror due to the stronger coupling be- 
tween the movable mirror and the cavity field. 

We thank Dr. M. Aspelmeyer for bringing the experi- 
mental paper [3^ to our notice and for interesting corre- 
spondence. We gratefully acknowledge support from the 
NSF Grant No. PHYS 0653494. 
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